RAPID STABILIZATION IN A SEMIGROUP FRAMEWORK 
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Abstract. We prove the well-posedness of a linear closed-loop system with 
an explicit (already known) feedback leading to arbitrarily large decay rates. 
We define a mild solution of the closed-loop problem using a dual equation 
and we prove that the original operator perturbed by the feedback is (up to 
the use of an extension) the infinitesimal generator of a strongly continuous 
group. We also give a justification to the exponential decay of the solutions. 
Our method is direct and avoids the use of optimal control theory. 



1. Introduction 

We consider a physical system which state x satisfies the Cauchy problem 

Bu(t), 




where A is a linear differential operator that models the dynamics of the system 
and B is a control operator that allows us to act on the system through a control u. 

The stabilization problem consists in finding a feedback operator F such that the 
solutions of the closed-loop problem 

x = (A + BF)x 

tend to zero as t tends to +00. 

For finite-dimensional systems, D. L. Lukes [18] and D. L. Kleinman [10] (see 
also the book of D. L. Russell [20]. pp. 112-117]) gave a systematic stabilization 
method thanks to an explicit feedback constructed with the controllability Gramian 

A := / e- tA BB*e~ tA * ' dt. 



The above matrix is positive-definite provided that (A, B) is exactly controllable 
(equivalently (— A*,B*) is observable). In this case, the feedback 

F := -B*^ 1 

stabilizes the system. 

Later, adding a suitable weight-function inside the Gramian operator, M. Slem- 
rod |21| adapted and improved this result to the case of infinite-dimensional systems 
with bounded control operators. More precisely, his feedback depends on a tuning 
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parameter u) > that ensures a prescribed exponential decay rate of the solutions. 
The weighted Gramian 

r T 

A w := / e - 2ujt e- tA BB*e- tA dt 
Jo 

is positive-definite if (A, B) is exactly controllable in time T (or (—A* , B*) is exactly 
observable in time T) . Then the solutions of the closed- loop problem provided with 
the feedback 

F := -5* Aj 1 

decrease to zero with an exponential decay rate being at least w i.e. there is a 
positive constant c such that 

\\x{t)\\ < ce- ut \\x \\, t>0, 

for all initial data xq, where || • |j denotes a norm on the state space. 

The problems that we have in mind are linear time-reversible partial differential 
equations (waves, plates...) with boundary control. These are infinite-dimensional 
problems and controlling only at (a part of) the boundary of the domain imposes 
that the control operator B is unbounded. This leads to difficulties in choosing the 
right functional spaces and the right notion of solution to have well-posed open-loop 
and closed-loop problems. 

J.-L. Lions |T7] gave an answer to the stabilization of such systems. His proof, 
using the theory of optimal control, is non-constructive and does not give any 
information on the decay rate of the solutions. By using a slightly different weight 
function in the above operator A u , V. Komornik gave an explicit feedback 
leading to arbitrarily large decay rates. His approach does not use the theory of 
optimal control : an advantage is that one does not have to use strong existence and 
uniqueness results for infinite dimensional Riccati equations. In fact, the weight 
function is chosen in such a way that A u is the solution of an algebraic Riccati 
equation. Formally, 

Aku + KA* + A U C*CA U - BB* = 0, 

where a definition of the operator C and the rigourous meaning of this equation 
will be given later. For a presentation of this method of stabilization, see also the 
books of V. Komornik and P. Loreti [13l pp. 23-31] (where a generalization of this 
method to partial stabilization is also given) and J.-M. Coron pp. 347-351]. 

Applications of this method to the boundary stabilization of the wave equation 
and the plates equation are given in [TT|. This method can also be used to stabilize 
Maxwell equations [112] and elastodynamic systems [T|. Moreover, numerical and 
mechanical experiments ([4], [5], [22]) have proved the efficiency of this feedback. 

In this paper, after recalling the construction of V. Komornik's feedback law and 
some results about the well-posedness of the open-loop problem (section[2]), we give 
a proof of two points that were not justified in |11| . 

• The first point (section [3]) is the well-posedness of the closed- loop problem 
with the explicit feedback introduced in |11| . Using the Riccati equation 
satisfied by A w , we introduce a "dual" closed- loop problem, which is easier 
to deal with because it does not involve the unbounded control operator 
B. Then we give a definition of the mild solution of the initial closed-loop 
problem and we prove in Theroem 13.11 that this solution satisfies a variation 
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of constants formula. To derive this formula, we adapt a representation 
formula of F. Flandoli [9] to the case of an algebraic Riccati equation. In 
Theorem 13.31 we prove that using a suitable extension A of A, the operator 
A — BB*A~ 1 is the infinitesimal generator of a strongly continuous group 
on the original state space. 
• The second point (section 2} consists in the justification of a formula, con- 
tained in Proposition 221 that is used in [11] to prove the exponential decay 
of the solutions. We recall at the end of the paper how this formula is used 
to obtain the exponential decay. 

The author thanks Professors I. Lasiecka and R. Triggiani for fruitful conver- 
sations on the subject of this work. 



2. A SHORT REVIEW OF THE CONSTRUCTION OF THE FEEDBACK AND OF THE 

OPEN-LOOP PROBLEM 

2.1. Hypotheses and notations. The state space H and the control space U are 
Hilbert spaces. We denote by H' and U' their duals and by 

J : U' — >■ U the canonical isomorphism between U' and U; 
J : H H' the canonical isomorphism between H and H' . 
Moreover we make the following hypotheses : 

• (HI) The operator A : D(A) C H — >• H is the infinitesimal generator of a 
strongly continuous group e tA on H. Q 

• (H2) B G L(U, D(A*)'), where D(A*)' denotes the dual0 space of D{A*). 
Identifying D(A*)" with D(A*), we denote by B* G L(D(A*),U') the ad- 
joint of B. This implies the existence of a number A G C and a bounded 
operator E G L(U,H) such that 

B* =E*(A + \I)*. 

• (H3) Given T > 0, there exists a constant C\(T) > such that 

f \\B*e- tA 'x\\l,dt <c x {T)\\xf H , 
Jo 

for all x G D(A*). In the examples, this inequality represents a trace 
regularity result (see |14]). It is usually called the direct inequality. 

• (H4) There exists a number T > and a constant c^iT) > such that 

C2(T)\\x\\ 2 a , < [ T \\B*e- tA *xf ur dt 
Jo 

for all x G D{A*). It is usually called the inverse or observability inequality. 



1 Thus his adjoint A* : D(A*) C H' — > H' is also the infinitesimal generator of a strongly 
continous group e tA = (e tA )* on H' . 
2 Provided with the norm 



||2 

Ud(A') 



D(A*) is a Hilbert space. Moreover, 

D(A*)CH' ffcfl(A')'. 
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Remark. Thanks to the assumptions (H1)-(H2), if the direct inequality in (H3) is 
satisfied for one T > 0, then it is satisfied for all T > 0. Moreover, the estimation 
remains true (up to a change of the constant in the right member) if we integrate on 
(— T, T). Extending this inequality to all x £ H' by density, the map t n- B*e~ tA x 
can be seen as an element of Lf oc (R;U'). 

2.2. Construction of the feedback. 



The operator A^. We suppose that the hypotheses (H1)-(H4) hold true (the 
number T > giving the observability inequality in (H4)) and we recall the con- 
struction of the feedback exposed in [11] , by defining a modified, weighted Gramian. 
We fix a number w > 0, set 

and we introduce a weight function on the interval [0,T U ] : 

j e~ 2u3S si < s < T 
6a,(s) := ( 2ue- 2 » T (T ul - s) si T<s<T u . 

Thanks to (H3) and (H4), 

{Kx,y) H ,W ■= [ e LJ {s)(JB*e~ sA 'x : B*e~ sA 'y) u ^ u ,ds 
Jo 

defines a positive-definite self-adjoint operator G L(H' , H ). Hence A„ is invert- 
ible and we denote by Aj 1 e L(H,H') its inverse. 

Actually the weight function e w has been chosen in such a way that the operator 
A w is solution to an algebraic Riccati equation. We are going to derive this Riccati 
equation because it will play a key role in the analysis of the well-posedness of the 
closed-loop problem and the exponential decay of the solutions. 



An algebraic Riccati equation. Let x.y € D{{A*) 2 ). We compute the integral 



(1) 



d 
ds 



eiLl {s)(JB*e~ sA "x,B*e 



D -sA*y\ 



U,U' 



ds 



in two different ways. Note that the quantity between the brackets is differentiable 
in the variable s thanks to the regularity of x and y, and the hypothesis (H2) made 
on B*. 

• On the one hand, as e w (T w ) = and e^(0) = 1, the above integral is 

-(JB*x,B*y) UtU ,. 

• On the other hand, by differentiating inside the integral, we obtain 

e' u (s){JB*e- sA x,B*e~ sA y)u,u>ds - / e w (s)( JB* e - sA A*x,B*e~ sA y)u,u>ds 



e LU {s)(JB*e- sA * x,B*e- sA ~ A*y) u , u ,ds 



-sA* 
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The formula 

(Lx,y) H :=- [ " el(s){JB*e- sA *AZ 1 x,B*e- l,A *AZ 1 y)u,u>ds 
Jo 

defines a positive-definite self- adjoint operator L g L(H) because 

Vs > 0, -e' u (s) > 2we w (s). 

We set 

C :=VIeL(H). 

For x,y £ H, we have 

(Lx,y) H = (Cx,Cy) H 

= (Cx,JCy) H ,H> 
= (x, C* JCy)H.H' 

where C* € L(H') is the adjoint of C. We can also remark the important relation 
between C and A J 1 : 

(2) C*JC > 2wAj x . 
Finally the second computation of the integral gives 

- (LAujX^^h - (A u A*x,y) H ,H> - (A LJ x,A*y) H . H > 
= - (CA^XiJCA^h.h, - (A u A*x,y) H ,H> - (A u x,A*y) H ,H'- 

Putting together the two computations, we obtain the following algebraic Riccati 
equation satisfied by A w : 

(3) (A UJ A*x, y) H .w + (A u x, A*y) H , H , 

+ (CA 0J x,JCA u y) HtH , - (JB*x,B*y) u , u , = 0, 

first for x,y E D((A*) 2 ) and then for x, y <E D(A*) by density of D((A*) 2 ) in D(A*) 
for the norm || • 

An integral form of the algebraic Riccati equation. We rewrite the Riccati equa- 
tion |3]) in an integral form, verified for x,y € H instead of x,y £ D(A*). Set 
x,y E D(A*). The equation ((3]) applied to e~ sA "x, e~ sA * y g D{A*) gives 

(4) {A^e-^'x, e- sA *y) H>w + (A w e-" A ' x, A*e' sA 'y) H M' 

+ (CA u e- sA 'x,JCA u e- sA 'y) H ,H' ~ (JB*e~ sA 'x, B*e~ sA ' y)u, w = 0. 
Integrating ^ bewteen and t gives the following integral form of the Riccati 
equation © : 



(5) (A u x,y) H> H' = (K>eT tA x,e~ tA y) H ,H' 

(CA LJ e- sA *x,JCA 0J e- sA *y) H M>ds+ / {JB* e - sA ' x,B*e~ sA ' y)u^d S . 
o Jo 

This relation remains true for x,y € H' by density of D(A*) in H' for the 

norm II • \\ H ,. 



^ This estimation is important for the proof of the exponential decay of the solutions. 
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Rapid stabilization. Now let us recall the main result of 

Theorem 2.1 (Komornik, |TTJ p. 1597]). Assume (H1)-(H4) for some T > 0. Fix 
lj > arbitrarily and set 

F := -JB*A- 1 . 

Then the operator A + BF generates a strongly continuous group in H and the 
solutions of the closed-loop problem 

x' = Ax + BFx, x(0) = x 

satisfy the estimate^ 

< \\xo\\ u e- ut 

for all xq G H and for all t > 0. 

2.3. Well-posedness of the open-loop problem. In this paragraph, we recall 
some results about the well-posedness of the open-loop problem 



(6) 



x'(t) = Ax(t) + Bu(t), te 

x(0) = .T , 



where u G L[ 2 oc (IR; U). We would like to define a mild solution of this problem that 
is continuous and takes its values in H . The dificulty comes from the fact that the 
control operator is unbounded and takes its values in the larger space D(A*)' . The 
next proposition will give an answer. 

Proposition 2.2 ([3, p. 259-260], [H p. 648] ). Fix T > and set 

z(t) : = / e^- s)A Eu{s)ds, —T < t <T. 

Then 



o 



• z(t) e D(A) for all -T<t< T; 

• \\{A + \I)z(t)\\H < k\\ u \\L 2 (~T,T;U) f or M —T < t < T, where k > is a 
constant independant of u; 

• (A + \I)z G C([-T,T];H). 

Definition. We define the mild solution of © as the application 

(7) x{t) = e tA x + (A + XI) f e^ t - s '> A Eu(s)ds 

Jo 

which is continuous on R with values in H . 

Remark. The relation ([7]) is a variation-of-constants-type formula. If B is bounded, 
this relation corresponds to 

x(t) = e tA x + [ e {t - s)A Bu{s)ds. 
Jo 



4 As it was already noted in pill , we have to consider this affirmation in a weaker sense. More 
precisely, we will see that this is true if we replace A by a suitable extension. 

° || ■ \\uj defined by \\x\\^ := (A^j 1 x,x) H i H , is a norm on H, equivalent to the usual norm 
thanks to the continuity and coercivity of • 

This result is due to I. Lasiecka and R. Triggiani who first proved it in the case of hyperbolic 
equations with Dirichlet boundary conditions (see |14| V 
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Moreover we can also write the relation by using the duality pairing : 

(x(t),y)u,H' = ( x o,e tA *y) H .H> + f (u{s),B*e^~^ A *y} UtU ,ds, 

Jo 

for all y € D(A*). 

We end this section by recalling a regularity result. It concerns the solutions of 
the open-loop problem in the dual space H' 



(8) 



y'(t) = -A*y(t)+g(t), te 
2/(0) = 2/q, 



where g € Aoc( R ; ■#')• Tnis time > tne 

source term does not involve any unbounded 
operator and the mild solution of ([SJ is defined by the "standard" variation of 
constants formula (see [HI p. 107] 



<<M y(t) = e- tA *y + f : e~^ A ' g(r)dr, 

Jo 

which is a continuous function from K to H'. Thanks to the direct inequality 
stated in (H3), we can apply the operator B* to the solution of the homogeneous 
problem associated to (|SJ) (put g = in @) and see this new function as an element 
of Lf oc (R; U'). Actually, this operation can be generalized to the solutions of the 
inhomogeneous problem (g can be ^ 0). We recall this result in the 

Proposition 2.3 ([9, pp. 92-93] , [16, p. 648] ). Fix T > 0. There exists a 
constant c > such that for all yo G D(A*) and all g 6 £ x (— T, T; D(A*)) we have 
the estimation 

\\B*y(t)\\l,dt < c(\\ yo \\ 2 H , + ||.g||i 1( _ T . T;C/0 ), 

-T 

where y is defined by ([9]) . By density, we can say that this estimation remains true 
for all initial data yo € H' and all source terms g e L X (-T,T;H'). 

3. WELL-POSEDNESS OF THE CLOSED-LOOP PROBLEM 

The aim of this section is to give a notion of solution to the closed-loop problem 

' x'(t) = Ax(t) - BJB*A~ 1 x(t), teK, 
x(0) = x . 



(10) 



As for the open-problem ((BJ), we can not use directly a variation of constants for- 
mula because of the unbounded perturbation (—BJB*A~ 1 ) of the infinitesimal 
generator A. 

Let us give the main idea for the well-posedness of (fTO)) . The Riccati equation 
§3§ can be rewritten formally as 

AA U + K U A* + k u C*JCK u - BJB* = 0. 

By multiplying the above equation on both side by Aj 1 , we get 

(11) K'^A + A* A- 1 + C*JC - A^ 1 BJB*A- 1 = 0. 



^ This result was firstly stated in 1141 in the case of hyperbolic equations with Dirichlet bound- 
ary conditions. 
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Now we multiply the operator A — BJB*A LJ 1 on the left by A^ 1 and on the right 
by A w to get 

AZ\A - BJS*A- 1 )A W = A^AA^ - A^BJB* = -A* - C*JCA Ml 
the last equality being a consequence of (fTT|) . 

Remark. The two operators A — BJB* and —A* — C*JCA UJ are (formally) conju- 
gated by the operator A w . 

The advantage of working with the conjugated operator is that the perturbation 
(— C* JCA U ) is bounded. We are going to analyze the well-posedness of the closed- 
loop problem (fTt)|) by using the solutions of the "conjugated" closed-loop problem 

JV (f) = -A*y(t) - C'JCKyfJt), i e K, 

( ] \v(0) = y o , 

whose well-posedness is already known. 

The perturbation being bounded, the operator — A*— C* JCA U , defined on D(A*) 
is the infinitesimal generator of a strongly continuous group V(t) on H' (see |191 p. 
22 and p. 76]). Moreover, for all t £ M. and all j/o G -ff' we have 

(13) V(t)y = e- tA 'y - / e ~^ A 'C* JCA u V(r)y dr. 

Jo 

Definition. Let Xq £ H. We define the mild solution of (fTU|) by 

C/(t)x := A^V^A^xq. 

Now we prove that this notion of solution is "coherent" with the closed-loop 
problem (JTUJ) in the sense that it satisfies a variation of constants formula, close to 
the one that we would formally use. 

Theorem 3.1. U(t) is a strongly continuous group in H whose generator is 

A v := A^(~A* - C* JCA^AJ 1 ; D(A V ) = A U D(A*). 
Moreover, it satisfies the variation of constants formula 

(14) (U(t)x 0) y) = (e tA x Q ,y) - [ (JB*Az 1 U(r)x ,B*e^ A "y)dr ) 

Jo 

for all xq £ H and y £ H' . 

Remark. The formula (fT4")l does not mean that A — BJB* Aj 1 is the infinitesimal 
generator of a group (or even a semigroup) but it justifies the choice of U (t) to define 
the mild solution of the closed- loop problem (JTUJ). To justify that the integral in 
is meaningful, see the remark after the Lemma just below. 

Remark. Formula (|14p can be rewritten as 

(15) U(t)x = e tA x -(A + XI) f e^ t ' r '> A EJB*AZ 1 U(r)x Q dr, 

Jo 

for all xq £ H. We can show (TTSl first for xq £ D(A*) and extend it by density. 
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The proof of Theorem 13.11 relies on the following representation formula of A u . 
Lemma 3.2. Set x,y £ H' and i 6 K. Then 

(16) (A u x,y} H ,H> = (A u V(t)x,e- tA *y) H ,H> + [ (JB*V{s)x, B*e~ sA * y)u,u>ds. 



Remark. The integral in the above formula is meaningful. Indeed the first part 
of the bracket defines an element of Lf oc (M.; U) because of (fTB")) and the extended 
regularity result stated in Proposition 12.31 The second part of the bracket defines 
an element of Lf oc (R; U') thanks to the direct inequality stated in (H3). 

Proof of Theorem \3.1\ At first, U(t) is a Co-group on H because it is the conjugate 
group (by A w ) of V(t) . The relation between the infinitesimal generator of V(t) 
and those of U(t) is also a general fact about conjugate semigroups (see [HI p. 43 
and p. 59]). 

To prove relation (|14p. we use relation in which we replace x by Aj 1 xo and y 
by e tA " y. Finally we use the definition of U(t), that is U(t)x = A w V{t)A~ x XQ. □ 



Proof of Lemma \3. 6 2[ F. Flandoli has proved in |9J a similar relation for the solution 
of a differential Riccati equation. We adapt his proof to the case of an algebraic 
Riccati equation. The proof contains two steps : at first, we use the integral form 
of the Riccati equation and the variation of constants formula for V (flUl) to 
prove relation (jTHJ) modulo a rest. Then we show that this rest vanishes. @ 

First step. Fix x,y £ H' and tel. From © and ((TS} we have 















~e~ tA 'x 


,e- tA *y)~ f (CA U 
1 Jo 


'e~ sA 'x 


, JCA UJ e~ sA "y)ds 


+ ./ 


f (JB* 




e~ sA 'x 


,B*e- sA 'y)d 


s 




(A u 


V(t) + 


/'-" 

Jo 


-r)A*C'JCKV(r)dr 


x,e- tA 'y) 



(CA W V(s)+ / e-t'-^'C'JCKVtfdr 
L Jo 

+ / (JB*\v(s)+ f e-( s - r)A 'C*JCA u V{r)dr 



x, JCA u e~ sA 'y)ds 



x,B*e' sA *y)ds 



={A w V{t)x,e- tA y)+ / ( JB*V(s)x, B*e~ sA y)ds + R. 
Jo 



Second step. To obtain relation (HHJ), we have to show that the rest R vanishes. 
To lighten the writing, we set 

g (r) := C*JCA w y(r)x £ C{R;H'). 



In order to simplify the notations, we will omit the name of the spaces under the duality 
brackets in this proof. 
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Let us rewrite the rest : 



R=(K / e-^ A 'g(r)dr,e- tA 'y) 



(CA u V(s)x, JCA 0J e~ sA 'y)ds 



(CA U 



,-{s-r)A' 



+ / (JB* / e -( s - r ^* 
JO Jo 
= :R\ — i?2 — ^3 + Ri- 

• We can also write R\ as 



g(r)dr,JCA^e" sA y)ds 



g(r)dr,B*e- sA y)ds 



Ri 



(A u e- {t - r)A g(r),e- {t - r > A e~ rA y)dr. 



The integrand of the above integral corresponds to the first term in the right mem- 
ber of by replacing x by C* JCA w V r (r)x = g(r), y by e~ rA y and t by t — r. 
Hence 

Ri= I (A u g(r),e- rA *y)dr 



+ 



(CA w e" sA g(r), JCA w e" sA e~ rA y)ds 



dr 



JB*e- sA * g(r), B*e~ sA r e~ rA y)ds 



A" . 



dr 



='-R\ + R' 2 - Ra- 

• We have 

R'l = i?2- 

The change of variable a := s + r and Fubini's theorem give 



R' 2 = 



(CA U 



g{r),JCA u e~ aA 'y)dadr 



Jr 
t pa 



(CA u e- ( -' T - r)A "g(r),JCA u e-' rA "y)drda 



>o Jo 
=R S . 



• II remains to show that R' 3 — R±. Difficulties arise since the operator B* 
is unbounded. The idea is to construct two approximations R' 3 (n) and R^n) for 
i?3 and i?4. We show that R 3 (n) = Ri{n) and that R' 3 (n) and R^(n) converge 
respectively to R' 3 and R4. 

Remark. A* is the infinitesimal generator of a Co-group in H' . Hence for sufficiently 
large n £ N, n lies in the resolvant set of A*. We set 

I n := n{nl- A*)- 1 £ L{H'). 

Then for all x £ H' , I n x £ D(A*) and I n x — > x as n — > 00 (see |19[ Lemma 3.2. 
p. 9]). Moreover, the sequence ||J n || is bounded from above independently of n. 
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Indeed, as A* is the generator of a group, it results from Hille-Yosida theorem ( |19[ 
Theorem 6.3 p. 23]) that for sufficiently large n € N, 

\\I n \\ = \\n{nI-A*)- 1 \\< "" 



n-P' 



where a and (3 are two positive constants. 
• For n sufficiently large, we set 



R' 3 (n) 



(JB*e~ sA ' I n g(r),B*e- sA e~ rA y)ds 



dr. 



The application between the duality bracket is measurable on the product space 
(0,t) x (0,i)- Moreover 



o Jo 

t r t 



JB*e- sA 'l n g(r),B*e 



V) 



dsdr 



< 



[ f (JB*e- sA 'l n g(r),B*e- 
Jo Jo 



y) 



dsdr 



(JB*e~ sA " I n g{r),B*e- sA * e~ rA ' y) 



ds 



dr 



(Fubini-Tonelli) 

— c / lls( r )ll-tf' ll e ~ r ' 4 Wii'dr (Cauchy-Schwarz, direct inequality) 

Jo 
<oo. 

Hence we can invert the order of the integrals in R' 3 (n). We get (first by doing the 
change of variable a := s + r) : 

rt ft 

R\ 



? sW=/ / (JB*e-( a - r ^ A * I n g(r),B*e- aA 'y)dadr 

JO Jr 

= 11 (JB*e-( a -^ A * I n g(r),B*e- aA ' y)drda. 
Jo Jo 

Finally, R' 3 {n) = J Q tp n (r)dr et R' 3 — J* tp(r)dr with the evident notations. For all 
< r < t, we have 



\fn(r) - v(r)| 



(JB*e 



* „ — s A 



[I n g(r)-g(r)},B*e 



y)ds 



< 



[JB*e 



* ^, — sA* 



[I n g(r)~g(r)],B*e 



v) 



ds 



<c\\I n g(r) — <?(?") || i?' || e r v\\h' (Cauchy-Schwarz and direct inequality). 

Hence (p n (r) — > <p(r) as n — > oo. Thanks to Cauchy-Schwarz, the direct inequality 
and because ||/n|| is bounded from above, we have 

\fn(r)\ < c\\I n g{r)\\ H ,\\e- rA *y\\ H , < c'\\g(r)\\ H ,\\e- rA *y\\ H ,. 
We can apply the dominated convergence theorem : R 3 (n) — > R 3 . 



The right side is measurable because it is the composition of two measurable functions, 
(we recall that B*e~ tA is well-defined in Lf (K; U')). In the left side we can replace B* by 
:= E*(A*^ + XI) where At, G L{H') is the Yosida approximation of A* (see [19]). For all 
x g D(A*), B^x — > B*x as k — > oo and B^ G L(H' , U'). Hence, the left-hand side of the duality 
bracket is measurable as a simple limit of continuous (hence measurable) functions on (0, i) X (0, t). 
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For sufficiently large n, we set 



Ri{n):= / {JB* / e^—*^' I n g(r)dr, B*e~ sA ' y)ds. 



But I n et e ( s r ^ A* commute and 

B*I n = E* (A + \I)*n(nI - A*)- 1 

= -nE* + (n 2 + n\)E*{nI - A*)- 1 e L{H'). 
Hence (see [2j p. 139] for interverting B* and the integral) 

B* f e-^ A * I n g{r)dr = f B* I n e-^ A * g(r)dr 
Jo Jo 

and 

R i (n)= [ [ (JB*I n e- {s - r)A *g(r),B*e- sA *y)drds = Ri } (n). 



Finally, for all < r < t, I n g{r) — > g(r) and ||/ n ff( r )|| < c ll<?( r )l|j the right member 
being integrable on (0, t). Thanks to the dominated convergence theorem, I n g — > g 
in L 1 (0,t;i?'). The estimation of proposition 12.31 gives 

B* f e - {s - r)A ' I n g(r)dr -> B* f er {s - r)A ' g(r)dr 
Jo Jo 

in L 2 (0, t; U 1 ). Hence R^{n) — > R4 and by unicity of the limit, R' 3 = R4. □ 

We can do a little better and link the infinitesimal generator of U(t) to the 
original operators involved in the closed- loop problem (|10[) . 



Theorem 3.3. The operator A : D(A) C H — > H can be extended to a bounded 
operator A : H — > D(A*). The operator 

A- BJB*A- 1 

coincides with Ajj (the generator ofU(t)) on D(Ajj) — K U} D{A*) i.e. 

Vx G D(Au), (A - BJB*A~ 1 )x Q = A v x e H. 

We first recall a classical extension result for the unbounded operator A to a 
bounded operator on H with values on the extrapolation space D(A*)' (see [HI pp. 
6-7] and pp. 21-22]). 

Lemma 3.4. The operator A : D{A) C H — » H admits a unique extension to an 
operator A G L(H, D(A*)'). Moreover this extension satisfies the relation 

(17) {Ax,y) D(A , y ,D(A') = (x,A*y) H>H ,. 

for all x £ H and y £ D(A*). 



Proof of Lemma \3.4\ The unicity of such an extension is the consequence of the 
density of D(A) in H. Provided with the norm || • ||_d(a*)j D(A*) is a Hilbert space 
and A* G L(D(A*), H'). We denote by A the (Banach-)adjoint of A* seen as a 
bounded operator between the Banach spaces D(A*) and H. Hence the definition 
of the (Banach-)adjoint gives 

A € L(H, D(A*)') 
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and for all x £ H and y £ D(A*), 

(Ax,y) D ( A ,y tD ( A .) = (x,A*y) H ,H' 

i.e. relation (fT"7l) is true. Moreover this new operator A defines extension of A i.e. 
the two operators coincides on D(A). Indeed from the above relation specialized 
to x £ D(A) C H, we get 

Vy£D(A*), (Ax,y) D(A *y tD(A *) = {Ax,y) H ,H> => Ax = Ax £ H. □ 

Proof of Theorem \3.3[ Instead of returning to the Riccati equation ([3]) , we are going 
to differentiate the variation of constants formula (fH)) . We know that for xo £ 
A UJ D(A*) = D(Au), the map 

1 1 y U(t)xo 

is differentiable and 

jU{t)xv = AuU(t)x . 

In particular if y £ H ' , then 

(^-U(t)x ,y) = (Ay (7(i)a;o, 2/). 

Differentiating (jT5J) with respect to t, we want to link the generator Ajj and the 
operator A — BJB*A~ 1 (a priori with values in D(A*)'). We remark that defining 
the domain of the latter operator is not clear. Let xq £ A UJ D(A*) and y £ D((A*) 2 ). 
First step. The map 

r ^ B*A- 1 U(r)x 

is continuous from K to U'. Indeed, setting yo :— A^xq £ D(A*), we have 
B*Az 1 U(r)x = B*AZ 1 A u V(r)A- 1 x 
= B*V(r)y 
= E*(A* +XI)V(r)y 

= E*{A* + C*JCA U - C*JCA U + \I)V(r)y 

= -E*(-A* - C*JCA u )V{r)y + E*(-C*JCA U + XI)V(r)y 

= -E*V{r){-A* ~ C*JCA u )y Q + E*(~C*JCA U + XI)V(r)y , 



the latter expression being continuous in r. 

Remark. On D(A*) = D(—A* — C*JCA U ), the operators V(r) et -A* - C*JCA U 
(generator of V(r)) commute (this is a general fact about semigroups) but a priori 
V(r) and A* do not commute. 

Second step. The map 

s^B*e sA 'y 

is differentiable on R with values in V . Indeed, as y £ D((A*) 2 ). we have (A* + 
XI)y £ D(A*) and 

B*e sA 'y = E*e sA '{A* + XI)y. 
The latter expression is differentiable with respect to s and its derivative is B*e sA A*y. 

!0 Again, when the name of spaces under the duality brackets are unnecessary, we omit them. 
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Third step. We deduce from the two previous steps that the map 

ti-> / (JB*AZ 1 U(r)x , B*e^ r)A ' y)dr 
Jo 

is differentiable on K and its derivative is the map 

t ^ f\jB*AZ 1 U(r)x , B*e (t - r *> A ' A*y)dr + ( JB*AZ 1 U(t)x , B*y). 
Jo 

It results that given two (regular) data xo G A bJ D(A*) and y G D((A*) 2 ), we can 
differentiate (U(t)xo,y) with respect to t and get 

j t (U(t)x ,y) =(e At x ,A*y)- J (JB*AZ 1 U(r)x ,B*e^ A "A*y)dr 
- (JB*AZ 1 U(t)x ,B*y). 
Replacing y by A*y in (| 14[) and reinjecting in the above relation, we obtain 

(18) j t (U(t)x ,y) = (U(t)x ,A*y) - (JB*AZ 1 U(t)x Q ,B*y). 

With the same regularity as above for xq et y, we have 

^{U(t)x 0) y) H ,w = (AuU(t)x ,y) H ,H> = (AuU(t)x ,y) D (A*y ,d(a*), 

where Ajj is the infinitesimal generator of U(t). We recall from Lemma [3.41 that 
A admits a unique extension to an aperator A E L(H, D(A*)'). Thanks to this 
extension we can link Ajj and A — BJB*A~ 1 . From (fT8|) and (fl~7|) we have, for 
xo G A U D(A*) and y G D{{A*) 2 ), 

j t (U(t)x ,y) H M> - (U(t)x ,A*y} H ,H> - (JB*Az 1 U(t)x ,B*y) u , u , 

= (AU(t)x ,y) D ^A*)',D{A') - {BJB*A~ 1 U{t)x Q ,y) D(A , ) ^ D(A « ) 
= ((A-BJB*AZ 1 )U(t)x ,y) D{ A*y,D(A*) 

= (AjjU(t)xQ, y) D(A')' ,D(A*) 

In particular, the latter equality is true for t = 0. Hence, given a fixed xo G 
A UJ D(A*), we have 

({A- BJB*A^ 1 )x Q ,y) D{A , ) ^ D(A , ) = (Aux ,y) D (a*)>,d(A'), 

for all y G D((A*) 2 ). This relation remains true for all y G D{A*) by density of 
D((A*) 2 ) in £>(A*) (for the norm || • \\ D ( A *))- Finally, 

V.t G A U D{A*) = D(Au), (A - BJB*A~ 1 )x a = A v x Q e H. □ 

Remark. With an unbounded control operator (i.e. B G L(U, D{A*)')), one can 
prove, through examples, that the domain of Ajj is not always included in the 
domain of A. Thus, in general, the extension A is necessary in order to link Ajj 
and A on D(Ajj) (i.e. we cannot omit the "tilde" in the above relation). This 
phenomenon does not appear with a bounded control operator (i.e. B G L(U,H)): 
in that case, we can prove that the spaces D{A) and D(Ajj) coincide. 
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4. Derivation of a representation formula for A^ 1 

AND EXPONENTIAL DECAY 

In this section, we give a justification to a representation formula for Aj 1 involv- 
ing the group U(t). This corresponds to the formula (3.11) in [TT]. We recall it as 
it is written in [11] : for all s,t € R, 

A-^Uit-syA^Uit-s) 
+ { U(t - s)*{C*JC + Az 1 BJB*Az 1 )U(t - s)dr. 

J s 

This formula is used in to prove the exponential decay of the solutions of the 
closed-loop system. Again, F. Flandoli derived an analog formula in the case of 
differential Riccati equations in [9]. We adapt his proof to the case of algebraic 
Riccati equations. 

We first prove a similar representation formula for A w . 

Proposition 4.1. For all x,y £ H' and t E K 

(19) (Aux, y) H ,H> = (A u V{t)x, V(t)y) H ,w 

+ f (JB*V(s)x,B*V(s)y) UtU ,ds+ f {CA u V{s)x,JCA u V{s)y) H , w ds. 
Jo Jo 

Proof. It relies on the representation formula (|16[) for A^ that we have already 
proved : for x, y G H', 

(A ul x,y) = (A L ,V(t)x,{e- tA *y})+ f (JB*V(s)x, B*{ e - sA * y])ds. 

Jo 

In the right member of the above relation, we replace e~ tA y and e~ sA y by using 
the variation of constants formula ([13)1 for V : 

(A u x,y) ={KV{t)x,V{t)y) 

+ {A u V(t)x, [ e^-^'^JCAuV^yds) 
Jo 

+ [ (JB*V(s)x,B*V(s)y)ds 
Jo 

+ f (JB*V(s)x,B* f e - {s ~ r)A ' 'C*JCA u V(r)ydr)ds 
Jo Jo 

= : r 1 + r 2 + T 3 + r 4 . 

But 

T 2 := f (A u V(t-s)V(s)x,e-^ A *C*JCA u V(s)y}ds. 
Jo 
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Thanks to (fl6l) . applied to V(s)x instead of x, C*CA ul V(s)y instead of y and t — s 
instead of t, we have 

T 2 = I (A UJ V(s)x,C*JCA ul V(s)y)ds 



f (A u V(s)x,C*JCA ul V(s)y)ds 
Jo 

/ / (JB*V(r)V{s)x,B*e~ rA 'C*JCA LJ V(s)y)drds 
Jo Jo 



10 JO 

= f {CA u V(s)x,JCA u V(s)y)ds 
Jo 

rt rt-s 



PZ pt — 8 

- (JB*V(r + s)x, B*e- rA ' C*JCA ul V(s)y)drds. 

Jo Jo 

The change of variable a := r + s in the last term gives 



T 2 = / (CA u V(s)x, JCA U] V{s)y)ds 
Jo 

-If (JB*V(a)x,B*e-^-^ A *C*JCA UJ V(s)y)dads 

Jo J s 

= f (CA u V(s)x, JCA u V(s)y)ds 
Jo 

- f f (JB*V(<r)x,B*e- ( -< y - s)A *C*JCA UJ V(s)y)dsd<T. 
Jo Jo 

Hence we have proved that 

(Kx,y) =(A u V(t)x,V(t)y) 

+ f (JB*V(s)x,B*V(s)y)ds 
Jo 

+ f (CA LJ V(s)x,JCA LJ V(s)y)ds 

(JB*V(s)x,B* f e-^-^ A 'C*JCA u V{r)ydr)ds 
Jo 

(JB*V(s)x,B*e- {s - r)A 'C*JCA u V(r)y)drds. 



t rs 



Jo 



We have already shown in the proof of Lemma 13.21 that the two last terms in the 
above relation cancel each other. Hence the relation is proved. □ 

Proposition 4.2. For all x, y E H and t EM. 

(20) (A" V y) H , tH = (A^Uftx, U{t)y) H ,. H 

+ f (JCU(s)x,CU(s)y) H ,, H ds + f (JB*A- x U(s)x, B* AZ 1 U{s)y) UtW ds. 



o Jo 
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Proof. We replace x by k u x and y by k u l y in the relation given by the 
Proposition 14.11 : 

(x,k-\j) =(K bJ V{t)k- 1 x,V(t)k- 1 y)+ f (JB*V(s)AZ 1 x,B*V(s)AZ 1 y)ds 

Jo 

+ f {JCKV{ 3 )kZ 1 x,Ck u] V(s)k- 1 y)ds. 
Jo 

Then, by definition of U , 

(A^ar.y) = {k- x U{t)x,U{t)y) 
+ I (JB*kZ 1 U{s)x,B*kZ 1 U(s)y)ds + [ (JCU(s)x,CU{s)y)ds. □ 



Remark. A simple change of variable implies that for all s, t G K and all x,y G H, 
(21) (kz'x, y) = (k- l U{t - s)x, U(t - s)y) 

+ { (JB*k- 1 U(T-s)x,B*k- 1 U(T-s)y)dT+ [ (JCU(t - s)x,CU(t - s)y)dr. 



Finally, let us recall the outline of the proof of the exponential decay of the 
solutions of the closeddoop problem (jTDJ). We denote by x{t) the mild solution of 
(HHD i-e. 

x(t) = U(t)x . 

Using the relation pi) with x = y = U(s)xq = x(s), we have 

(AJ 1 ^),^)) =(kZ 1 x(t),x(t)) 

+ { (JB*k- 1 x(T),B*k- 1 x(T))dT + { (JCx{T) 7 Cx(T))dT. 



Let < s < t. The estimation @ between C and A^ 1 and the positiveness of the 
second term of the right member in the above relation yield 

(Aj 1 a:( S ),x( S )) > {k- 1 x{t),x{t))+2u> f (k^x^), x(j))dr. 



A Gronwall-type lemma (see [TT| p. 1599]) gives 

< Mte- 2 ^ Vt>0. 
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